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Abstract: We use several diverse parameterizations of diffractive parton distributions, 
extracted in leading twist QCD analyses of the HERA diffractive deep inelastic scattering 
(DIS) data, to make predictions for leading twist nuclear shadowing of nuclear quark and 
gluon distributions in DIS on nuclei. We find that the HERA diffractive data are sufficiently 
precise to allow us to predict large nuclear shadowing for gluons and quarks, unambiguously. 
We performed detailed studies of nuclear shadowing for up and charm sea quarks and gluons 
within several scenarios of shadowing and diffractive slopes, as well as at central impact 
parameters. We compare these leading twist results with those obtained from the eikonal 
approach to nuclear shadowing (which is based on a very different space-time picture) 
and observe sharply contrasting predictions for the size and Q 2 -dependence of nuclear 
shadowing. The most striking differences arise for the interaction of small dipoles with 
nuclei, in particular for the longitudinal structure function F^- 
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1. Introduction 

The high-energy scattering of leptons and hadrons from nuclear targets offers a unique pos- 
sibility to study the distribution of quarks and gluons in nuclei, i.e., large parton densities 
in a large volume. It is now firmly established by several experiments (EMC and NMC 
at CERN, a number of experiments at SLAC, E665 at Fermilab) that at small values of 
Bjorken x, 5 x 10~ 3 < x < 0.03 — 0.07 (for Q 2 > 1 GeV 2 ), the inclusive nuclear structure 
function, F^, for DIS on a nucleus with A nucleons, is smaller than the incoherent sum 
of the nucleon structure functions, AF^ ■ This phenomenon is called nuclear shadowing 
(for recent reviews see |l]]). It is important to emphasize that the available data, which 
admittedly has a rather limited kinematic range, Q 2 <4 GeV 2 , demonstrate no significant 
Q 2 -dependence of nuclear shadowing. In addition, nuclear shadowing of the antiquark dis- 
tributions has been observed in the Drell-Yan process (2|, in the range 0.02 < x < 0.05 
and 16 < Q 2 < 35 GeV 2 . In other words, the data are consistent with the dominance of 
leading twist shadowing. 

Why is it still interesting to work on the subject of nuclear shadowing in DIS on nuclei 
after the topic has been investigated so thoroughly? There are several reasons. 

Firstly, an accurate evaluation of nuclear shadowing is vitally important for both light 
and heavy nuclei to pin down the nuclear parton distributions. These form the boundary 
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conditions for hard processes in electron-nucleus DIS experiments proposed in the USA (the 
electron-ion collider (EIC)) || and at DESY 0], in nucleus- nucleus and nucleon- nucleus 
collisions presently under way at BNL, and in proton- nucleus and nucleus- nucleus collisions 
planned for the LHC at CERN. 

Secondly, it is important to understand the intricate interplay of perturbative and 
non-perturbative effects in nuclear shadowing and its role in the parton distributions of 
nuclei at small x. Understanding this interplay is a key element in determining when the 
leading twist approximation (implicit in the DGLAP evolution equations) breaks down for 
sufficiently small x. 

Thirdly, because of the deep connection between nuclear shadowing and diffraction in 
DIS on a nucleon established by V. Gribov ||, accurate measurements of nuclear shadowing 
should help to discriminate between competing models of diffraction. 

The dynamics of nuclear shadowing in DIS at high energies is most transparent in the 
target rest frame where the virtual photon-nucleus (nucleon) scattering, 7* (g) + A{Pa) — ► 
X, is a time-ordered three-stage process. 

Firstly, the photon may be considered to fluctuate into a linear superposition of 
hadronic components labelled by the basis set \h k ). The components consist of quarks, 
antiquarks and gluons at small relative transverse distances (and their associated color 
fields) and of hadronic bound states at large transverse distances. For large center-of-mass 
energy, W (W 2 = (q + Pa) 2 ), i-e., small Bjorken x AQ 2 /W 2 , the fluctuation occurs a 
long distance upstream of the nuclear target. The incoming virtual photon (more precisely, 
the hadronic part of the virtual photon) state I7*) may be expressed generically as follows 

| 7 *> = I>fclM, (i-i) 
k 

in which the index "fc" should be understood to be a label which fully specifies the hadronic 
state (i.e., it includes information about momentum fractions and transverse momentum 
of the constituent partons, their helicities, etc.) and c k = (h k \j*). 

Secondly, the hadronic configurations \h k ) interact strongly with the nucleus 

= y,4'frj£*:' c fc , (1.2) 

where the subprocess cross section, d"^ k i, is related via the optical theorem to the imaginary 
part of the forward scattering amplitude for the process h k + A — > h k > + A normalized such 
that a kk = &h k A, the total /ifc-nucleus cross section. 
Finally, the hadronic final state X is formed. 

In general, each state h k suffers attenuation: at sufficiently high energies the states 
\h k ) interact coherently with several nucleons of the target. These multiple interactions 
decrease a kk , relative to the sum of a kk , of the individual nucleons, which can be translated 
to the corresponding inclusive structure functions as F^/iAF^ 1 ) < 1. Thus, the coherent 
interaction of the photon with several nucleons of the target nucleus invariably results in 
nuclear shadowing. 
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Equations (|L1,1.2) are very general. Any particular realization requires some mod- 
elling of the hadronic content of the (virtual) photon and of the interaction with the target. 
There are several models in the literature including Generalized Vector Meson Dominance 
(GVMD) models ||, two-phase models including vector mesons and the Pomeron |7j or 
other scaling contributions [||, and the QCD-improved aligned jet model In this paper, 
we review and compare the two most recent approaches to nuclear shadowing: the leading 
twist and eikonal approaches. 

The leading twist approach was developed in a recent work by Frankfurt and Strikman 
|1C], who attempted to place the treatment of nuclear shadowing on firmer ground. The 
authors generalized the Gribov theory of nuclear shadowing in hadron-nucleus scattering 
H to individual quark and gluon nuclear parton densities in light nuclei, probed in DIS 
on nuclei. As in the Glauber model, nuclei were treated as a dilute nucleon gas. The 
connection between nuclear shadowing in inclusive DIS on the nucleus and diffraction in 
DIS on the proton, the factorization theorem for hard diffraction [|ll|] and high accuracy 
DIS diffractive data allowed the leading twist (LT) component of the nuclear shadowing 
correction to nuclear parton distributions to be evaluated reasonably well. One immediate 
and very distinct consequence of the theory is that, at some initial scale Qq (Qq = 2 GeV 
in |i"0|]), nuclear shadowing for gluons is much larger than that for quarks (by a factor of 
about two or three for the lightest nuclei) because the diffractive structure functions are 
dominated by gluons. 

Another popular approach to nuclear shadowing is the eikonal approximation. This 
method is based on the assumption that the interaction between the fluctuation \hk) and 
nucleons of the target leaves the fluctuation unchanged (i.e., the interaction is diagonal 
in the appropriate variables). This allows the subprocess to be iterated and summed 
so that the \hk) -nucleus total cross sections, Oh k Ai entering Eq. ( |1.2| ) via o^^/, can be 
expressed in an exponential form via the total (Jh k N cross sections. This procedure is 
termed eikonalization, hence the name of the approximation. 

The profound conceptual difference between the leading twist and eikonal approaches 
is readily revealed by including QCD evolution. Indeed, in QCD quark-antiquark pairs and 
gluons are radiated by the partons constituting the fluctuations \hk) before and during the 
interaction with the target, so that the fluctuations \hk) necessarily mix. This is included 
in the leading twist approach and ignored in the eikonal approximation. As a result, the 
two considered approximations predict a different size and Q 2 -dependence of the shadowing 
correction. For example, the longitudinal structure function is very sensitive to the 
mixing of the qq and qqg fluctuations of the incoming photon and thus can serve as a 
sensitive observable to distinguish between the leading twist and eikonal approximation 
predictions. 

At the outset, it may be useful to clarify the relation of this work to a number of works 
on the small- x scattering off nuclei (for nice pedagogical reviews and references to original 
papers, see e.g., [12|). We are working in a kinematic region in which the leading twist 
contribution dominates diffraction and one has to take into account logs of Q 2 and one or 
two logs of xq/x in the evolution of the parton distributions. This can be accounted for by 
the NLO DGLAP evolution equations. The approaches of 12] focus on a resummation of 
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the logs of xq/x (cf. the BFKL approximation) which is not necessary for this kinematic 
region and neglect certain logs of Q 2 . 

This paper is organised as follows. In section [2] we review the leading twist approxima- 
tion for nuclear shadowing and give predictions for the quark and gluon nuclear densities 
and for the nuclear structure functions F^, using as input a broad range of leading twist 
diffractive parton distribution functions (DPDFs), which successfully describe the HERA 
hard diffraction data. We also review the eikonal approximation approach to nuclear shad- 
owing and discuss various serious problems associated with it. We dedicate section ||to the 
longitudinal nuclear structure function, F£, which is predicted to be very different within 
the two approaches, and indicates general qualitative patterns of contrasting behaviour for 
hard processes. Discussions and conclusions are presented in section ^. In appendix [A] we 
specify and compare the different parameterizations of the DPDFs and observe that, in 
spite of significant differences in the details of these, they all lead to fairly similar prob- 
abilities for quark-induced and gluon-induced diffraction and nuclear shadowing at small 
x. We also observe that all of these parameterizations lead to a very large probability of 
leading twist diffraction in gluon-induced processes, which is comparable to the black-body 
limit of 50% for x < 10~ 4 at Q 2 ~ 4 GeV 2 . Finally, the parameterization of the nuclear 
one-body density for a range of nuclei considered is given in appendix M. 



2. The leading twist and eikonal approaches to nuclear shadowing 

In this section we review both the leading twist and eikonal approaches to nuclear shad- 
owing (subsections 2.1 and 2.2). A critical comparison of the approaches is presented in 
subsection 



2.1 The leading twist approach 

The basis of our present day understanding of nuclear shadowing in high-energy scattering 
on nuclei was formulated in the seminal work by V. Gribov |J]. The key observation 
was that, within the approximation that the radius of the strong interactions is much 
smaller than the average inter-nucleon distances in nuclei, there is a direct relationship 
between nuclear shadowing in the total hadron- nucleus cross section and the cross section 
for diffractive hadron- nucleon scattering. While the original derivation was presented for 
hadron-deuterium scattering, it can be straightforwardly generalized to lepton-nucleus DIS. 

2.1.1 The master formulae for leading twist nuclear shadowing 

In this subsection we aim to motivate a master formula which relates nuclear shadowing 
to diffraction off a single nucleon and includes the scattering off an arbitrary number of 
nucleons. The main result of the work by Frankfurt and Strikman is the observation 
that the Gribov theory can be generalized to calculate the leading twist component of 
nuclear shadowing for each nuclear parton distribution separately. We shall review the 
principal steps of the derivation below. The starting point is the generalization of Gribov's 
result to DIS of leptons on deuterium. The nuclear shadowing correction to the deuteron 
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structure function, 5F? = F% + F 2 " - F 2 D , where F 2 D , F 2 P , F 2 n are the deuteron, proton, 
neutron structure functions, respectively, can be presented in the form 

1 2 r 

5F?{ Xl Q 2 ) = J dk 2 dx 1P F^ 4) ((3,Q 2 ,x 1 p,k 2 )F D (4k 2 +Ax 2 ]P m 2 N ) . (2.1) 



Here F^^ is the diffractive structure function of the proton in terms of the usual diffractive 

l\ - k 2 ) « Q 2 /(Q 2 + M 2 X ) 



variables, Q 2 = -q 2 , = -q 2 /(2q ■ k) = Q 2 /(Q 2 + M\ — k 2 ) « Q 2 /(Q 2 + M\) and 



xp = 2q ■ k/(2q ■ p) = x/f3 « (Q 2 + M|-)/(2p ■ g); &t is the transverse momentum transfer 
to the nucleon; r/ is the ratio of the real to imaginary parts of the diffractive scattering 
amplitudes 1 for the reaction 7* + N — > X + N; Fz> is the deuteron electromagnetic form 
factor 2 . 

The relevant Feynman graph for the forward scattering amplitude of the double scat- 
tering diagram is shown in Fig. |l], which also serves to define the four-momentum flow. In 
particular, the momentum transfer to the proton k = p' — p is mostly transverse 3 , so that 
k 2 ~ —k 2 . In order to find the shadowing correction to the total virtual photon-deuteron 
cross section (and hence SF®) one needs to find the imaginary part of the forward ampli- 
tude of the diagram in Fig. |]. 

We can also obtain Eq. (2.1) in a somewhat different way than in the original Gribov 



paper by using the Abramovsky-Gribov-Kancheli cutting rules [14] . Instead of considering 
directly the imaginary part of the total rescattering amplitude, one can first consider the 
cross section corresponding to the diffractive cut of the amplitude in Fig. || (see Fig. 1 of 
flCfl ). In this case, the answer is given by the same diffractive cross section that enters 
lepton-nucleon diffractive scattering, i.e., the diagram is proportional to |ImA| 2 + |ReA| 2 , 
where A is the amplitude for diffractive lepton-proton scattering. One then compares 
the expression for the imaginary part of the forward amplitude of the double scattering 
diagram to that with the diffractive cut. Using the cutting rules, one observes that the 
former is proportional to |ImA| 2 — |ReA| 2 , while F®^ is proportional to [Im^ 2 + |Rej4| 2 . 
This explains how the answer for can be expressed in terms of the nucleon diffractive 
structure functions, -FrP^, as well as the presence of the factor (1-rf) / (l+rf) in Eq. (fjj). 

The deuteron form factor arises from the non-relativistic calculation of the overlapping 
integral of deuteron wavefunctions tpjj. Indeed, when one nucleon with a three-momentum p 
receives a three- momentum k=(kt, ks ~ itinxjp), the other nucleon with a three- momentum 
— p receives a three- momentum —k (see Fig. |j]). Hence, the lower (nuclear) part of the 
diagram in Fig. || is proportional to J d 3 p i^d{p) ip(p + k) = Fu(Ak 2 ) because the electro- 
magnetic form factor is defined as F£)(q 2 ) = J d?p iPd(p) ip(p + <?/2). When the argument 
of the form factor is large, i.e., when either xjp > l/(mjvf\D) or kt S> I/Vd, the damping 
of the form factor annihilates 5F® (ro is the average electromagnetic deuteron radius). 



lr The factor associated with a non-zero 77 is absent in the original deriv ation of since the intercept of 
the Pomeron trajectory ajp(O) was assumed to be equal to unity, see Eq. (2.3). 

2 For simplicity we give the expression for the spin- less deuteron. For the spin-one deuteron the relevant 
form factors were calculated in 113]. 

The transverse plane is defined relative to the virtual-photon nucleon center of mass frame. Given that 
the momentum is mostly transverse, to a very good approximation, we use t and A: 2 interchangeably. 
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Equation ([O]) can be generalized to DIS on any nucleus. The shadowing correction 
to the nuclear inclusive structure function, which arises due to scattering on any pair of 
nucleons, 5F^ , can be written in the form 



5R 



A(2) 



A(A-l) 



167^ 



(1 - irjf 
l + rj 2 



d 2 b 



dzi 



X P,0 



dxjp x 



F 2 D ^(P,Qi, X]P ,k?) 



p A (b,z 1 )p A (b,z 2 )e i ^ m ^-^ 



(2.2) 



Here pa is the nuclear one-body density per nucleon normalized such that J d 3 fpA(r) = 1; 
Z\ and Z2 are the longitudinal coordinates of the two nucleons involved in the scattering; b is 
the impact parameter of the projectile with respect to the center of the nucleus. The upper 
limit of integration, xjp^q will be discussed later in the text. A non- vanishing longitudinal 
momentum transfer to the target is taken into account by the factor 4 exp(ixjpm7v(<zi — £2))- 
The implied suppression for large x jpm jy is implicitly included in Eq. (|2.1|) by the Axjpmj^ 
term in the argument of the deuteron form factor. 

We determine the ratio of the real to imaginary parts of the diffractive scattering 
amplitude, 77, by assuming a universal exchange and using the analytical properties of 
the scattering amplitude as a function of energy, as suggested by Gribov and Migdal [16]. 
Applying the Gribov-Migdal expression to our case, we determine rj as a function of the 
intercept of the Pomeron trajectory, ajp(O), 



7T 



V 



7T dlnyl 



2<91nl/x 2 



(a.p(0) - 1) 



(2.3) 



4 This factor is analogous to the one which was derived within the vector meson dominance model of 
photoproduction and electroproduction (for a review see (l^]). In the context of DIS it was first introduced 
in the second reference of Q. 
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Different values for ajp(O) obtained by various groups correspond to different values for 
77. In our analysis, we used rj = 0.22 in conjunction with the ACWT parameterization of 
DPDFs jlTj (ajp(0) = 1.14) and rj = 0.32 in conjunction with the HI [18] parameterization 
(a F (0) = 1.20). 



The integration over kt is absent in Eq. (2.2) since it is assumed that the slope of the 
elementary diffractive amplitude for the reaction 7* + iV — > X + N is much smaller than the 
slope of the nuclear electromagnetic form factor 5 . This allows one to replace the diffractive 
nucleon structure function by its value at k\ = 0. 



Equations (2.1) and (2.2) include only the spin non-flip part of the diffractive cross 
section, which is proportional to F 2 D ^ . Spin-flip effects in diffraction are highly suppressed 
at small t by the structure of the Pomeron vertex and thus can be safely neglected. 

The crucial observation made in [lCj was that Eq. (|2.2[) can be generalized to nuclear 



parton densities. Indeed, the QCD factorization theorems for inclusive DIS [19] and hard 
diffraction JO]] allow one to express the inclusive and diffractive structure functions, en- 
tering the left and right hand sides of Eq. (|2.2; ), respectively, as a convolution of the same 
hard scattering coefficients with the corresponding, i.e., inclusive and diffractive, parton 
distributions. Equating the terms in front of each hard scattering coefficient, one imme- 
diately arrives at the key expression for the shadowing correction to the nuclear parton 
distributions, 5f^ A (x, Q 2 ), 



5ff^Q 2 ) = ^-^l^Re 



-j— j- 2 / d'h I di\ I </:-. / d.r n > < 



1 / az 2 

x <y z\ J x 



fj/ N (P,Q ,xjp,0) Pa(Mi) p A (b,z 2 ) e 



ix P m N (z 1 -Z2) 



(2.4) 



Here j indicates a generic parton label (i.e., a gluon, or a quark of a particular flavour); fjj N 
is the diffractive parton distribution function (DPDF) of the nucleon 6 , for parton of type j. 
Equation (2.4) is the main result of |l0| and demonstrates the intimate connection between 



the nuclear shadowing correction to nuclear inclusive parton distribution functions (PDFs) 
and the DPDFs of the nucleon. Since fjj N obeys the leading twist DGLAP evolution 

equation, the Q 2 -evolution of &f^\ is also governed by DGLAP, i.e., it is by definition 
a leading twist contribution. This explains why the approach of [10] can legitimately 
be called the leading twist approach. However, the validity of Eqs. (2.1,2.2,2.4) does 
not require the absence of the higher twist effects in diffraction. Hence, as soon as a 
particular parameterization of the diffractive cross sections fits the data, one can use this 
parameterization to calculate the shadowing correction to the structure functions, using 
Eqs. (2.1,2.2,2.4), without addressing the question of decomposition of the diffractive cross 
section over twists. 



5 This is analogous to the optical model approximation of the Glauber model which is known to work 
very well for A > 12. 

6 We do not distinguish between protons and neutrons since we consider small- a; scattering in the vacuum 
channel. 
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Ignoring for a moment the contribution to nuclear shadowing arising from the interac- 

(2) 

tions with three and more nucleons of the target, the connection of SfjA t° nuclear {fj/A) 
and nucleon (/j/jv) PDFs is given by 

fj/A = Af j/N - Sfj% , (2.5) 

and is depicted graphically in Figs. |2| and [3| for the up quark and gluon distributions, 
respectively. 

u A =A u N -A(A-l)/2 shadowing 





a) u A 



b) u N 



c) shadowing 



Figure 2: Feynman diagrams corresponding to Eq. ( p.5[ ) representing the nuclear quark distribution 
(a) as an incoherent sum of the quark distributions in the nucleons (b) minus the nuclear shadowing 
correction (c). Note that one needs to take the imaginary part of the forward amplitudes in order 
to find the structure functions. 



However, Eq. ( |2.4D describes only one piece of the shadowing correction, namely the 
part which arises from the interaction with any two nucleons. Since the strength of the 
interaction with several nucleons is large, for a sufficiently heavy nucleus the interaction 
with three and more nucleons plays an important role. Such an interaction cannot be 
expressed directly in terms of diffraction on the nucleon and has to be modelled. We 
shall use a scheme inspired to a certain extent by the quasi-eikonal approximation, which 
successfully describes soft total hadronic cross sections of high energy hadron-nucleus in- 
teractions. In this approximation, the amplitude for the forward virtual photon-nucleus 
scattering is given by graphs containing virtual photon- nucleon cross sections of two kinds. 
These are the cross sections of diffractive scattering and the production of inelastic inter- 
mediate states, and cross sections of "elastic" rescattering of the produced inelastic states. 



In the formalism of cross section eigenstates [2C], where one introduces the notion of a 
distribution of cross sections for the scattering state, the amplitude for the scattering off 
k nucleons is proportional to (c ) where (...) denotes averaging over the scattering state 
wavefunction. The quasi-eikonal approximation corresponds to the approximation 

(a*) 



a 



k-l 
eff 



(2.6) 
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g A =A g N -A(A-l)/2 shadowing 





a); 



b) g N 



c) shadowing 



Figure 3: Feynman diagrams corresponding to Eq. (^j) representing the nuclear gluon distribution 
(a) as an incoherent sum of the gluon distributions in the nucleons (b) minus the nuclear shadowing 
correction (c). Note that one needs to take the imaginary part of the forward amplitudes in order 
to find the structure functions. 



One can interpret a e ff as the average strength of the interaction of the produced config- 
urations. In this approximation the total cross sections can be "eikonalized" or summed 
to infinity in an exponential form in order to present the answer in a compact analytical 
expression (cf. |15|1 ). Introducing the full nuclear shadowing correction to fj/ A , ^fj/A: 
where $fj/ A = ^fj/N ~ fj/A> the quasi-eikonal approximation for dfj/ A reads 



Sf j/A (x,Q 2 ) = A{A „ l) l^Re 



(1 ~ irjf 



oc 



d z b / dz 



CO 



i / dz2 / dxp x 

— X J Z\ J X 



f? /N (P, Q 2 ,X F , 0) PA (b, Zl ) PA (b, Z2 yz P rn N ( Z1 - Z2)e -(A/2 ){ l- irlKff J* dzpA ( z) 



The rescattering cross section cj 3 e * , is discussed further in appendix 



• (2-7) 



Equation (2.7) implies that the rescattering cross section for the interaction with three 
and more nucleons, fg//' ^ s ^ ne same as the rescattering cross section for the interaction 
with two nucleons. In general, this is not true since hadronic fluctuations with very differ- 
ent cross sections contribute to the total cross section of the interaction of a hard probe 
(virtual photon, VK-boson, etc.) with a parton of the nucleus, and this broad dispersion 
over cross sections should be kept in mind. However, for nuclei with realistic finite A, 
numerical studies with models which include effects associated with the distribution over 
cross sections demonstrate that deviations from the quasi-eikonal approximation are small 



1 21, 22]. Similar conclusions were reached in [23] for the case of nuclear shadowing, where 
an even broader range of possible cross section fluctuations was considered. 

We used Eq. ( |2.7| ) in order to estimate Sfju at some initial Q 2 = Qq (Qq = 4 GeV 2 
in our case). The result was used as an initial condition for the QCD evolution of fj/ A 
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to higher scales Q 2 . We would like to stress that it would be incorrect to apply Eq. ( |2.7| ) 
at any arbitrarily large Q 2 : as Q 2 increases, the QCD radiation of gluons and qq pairs 
increases the dispersion of cross section fluctuations and hence makes the quasi-eikonal 
approximation less justified. 



2.1.2 The transition between shadowing and antishadowing: setting xjpo- 

The upper limit of integration over xjp in Eqs. (2.2, 2.4, 2.7), xp$i deserves a special discus- 



sion. In this subsection we motivate our choices for the numerical values of this parameter 
for the various parton species and explain that in our model this parameter effectively acts 
as the transition point from the shadowing region at low x to the antishadowing region 



at larger x. In [ 10 1 , it was assumed that xjp^ = 0.02. This choice of xjp^q was motivated 
by two observations. Firstly, only for xjp < 0.02 is it valid to assume that the Pomeron 
contribution dominates diffraction, i.e., that sub- leading Reggeon exchanges can be safely 
neglected. Secondly, larger xpp means that Eq. ( |2.4| ) can be applied to larger Bjorken x. 
However, for x > 0.02, nuclear parton distributions are expected to receive an additional 
contribution from nuclear antishadowing (the mechanism leading to the experimentally 
observed enhancement of nuclear inclusive structure functions), which is not included in 
Eq. ( |2,4| ). This makes the application of Eq. ( |2.4p for x > 0.02 not well motivated. 

The HERA diffractive data indicate that Pomeron exchange dominates diffraction for 
xjp < 0.01. At higher xp it becomes necessary to take into account sub-leading Reggeon 
exchanges, although the "Pomeron term" appears to be giving the dominant contribution 
for xjp < 0.03. At the same time the separation of the diffractive cross section into these two 
contributions in this kinematic region is definitely not unique and leads to some uncertain- 
ties. However, the sensitivity of the predictions resulting from Eq. ( |2.4| ) to the particular 
choice of xp$ is greatly reduced by the factor associated with the nuclear wavefunction, 
i.e., by the suppression factor paQ>, Zi)pA(b, z-i) exp(ixpmjsr(zx — 22))- 

Moreover, to perform a self consistent description of the parton densities for the whole 
range of x one needs to incorporate effects of enhancement of the parton densities at 
higher x. For the gluon such an enhancement follows from an energy-momentum sum rule 



analysis of the DIS nuclear data [24, 25] which revealed that the light-cone fraction carried 



by gluons in a nucleus and in a free nucleon is practically the same. For the valence quarks 



such an enhancement follows from the baryon charge sum rule [24]. For the antiquark/sea 
channel there are no general arguments to determine the effect of the enhancement and 
in fact no enhancement is observed in the Drell-Yan process 7 0. Hence we followed a 
somewhat simplified scenario of |m| , [29] which assumes that the enhancement is present 
only in the gluon and valence quark channels. The absence of the enhancement for the sea 
in principle allows one to take xp^ arbitrarily large. However, in our numerical analysis 
we set xjpfl = 0.1 for the sea. We explicitly checked that the answer is insensitive to the 
particular choice of xp t o due to the nuclear wavefunction suppression and the dominance 
of small- x values for the sea quark distributions. 

7 We note in passing that the observed pattern of the yl-dependence of nuclear parton densities at 
x ~ 0.1 — 0.2 may be related to the origin of nuclear forces |24, M. 
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We take xp 7 o = 0.03 for the gluons. For the gluon case we need to introduce a 
transition from the shadowing to the antishadowing regions for some particular value of 
x = xtr, which in general should depend on A. Arguments were given in p4| to the 
effect that the enhancement should be restricted to the region of x < 0.2, and could start 
somewhere in the range xt r £ [0.02,0.05]. In principle we could keep the same Xjp,o = 0.1 
as for the sea quarks and introduce an enhancement term for x > xt r to compensate for 
the shadowing term and to ensure an enhancement. As one can see from the structure 



of Eqs. (2^, 2A, ^j), the shadowing correction vanishes when x = xjp,o, and is strongly 
suppressed close to it. Hence we find that for our numerical analysis it is simpler to take 
xpfi = xtr and to neglect the enhancement term for x < xt r - 

The transition point between the shadowing and antishadowing regions for gluon dis- 
tributions in nuclei has never been measured experimentally. The best one can do is to 
infer some information about Xf r for the gluons indirectly, using QCD evolution. Our 



choice of xjp^ for the gluon distributions was motivated by the analysis in [26] of the NMC 
high statistics data on the ratio of structure functions F2 of Tin ( 118 Sn) and Carbon ( 12 C) 
1 27]. This analysis indicated that the transition between shadowing and antishadowing 



occurs for 0.023 < x < 0.035 and Q 2 4 GeV 2 . Thus, we use xp$ = 0.03 for the gluon 
distributions in our analysis. However, a choice of xjp,o does not fix entirely the shape of 
the gluon distribution ^(x). We assumed a quadratic polynomial form for the antishad- 
owing correction to qa{x) on the interval 0.03 < x < 0.2. The three free coefficients of the 
quadratic fit were chosen such that the antishadowing correction vanishes at x = 0.03 and 
x = 0.2 and that the momentum fraction carried by gluons bound in nuclei is the same 
as that in the free nucleon [25]. A similar shape of the resulting gluon enhancement was 



assumed in [25, 2 



It is important to mention that, in this work, we make predictions for the shadowing 
corrections to the gluon and sea quark distributions in nuclei and do not concern ourselves 
with the valence quarks. In our numerical analysis, we keep only the Pomeron exchange 
contribution to , which has vacuum quantum numbers. Thus, nuclear shadowing for 
the valence quarks in nuclei, which is associated with the Reggeon exchanges with non- 
vacuum quantum numbers, is not considered and effects of the valence quark enhancement 
for x ~ 0.1 are neglected. In general, a complete model for the leading twist parton 
distributions in nuclei would require conservation of the baryon number and momentum 
sum rules || in order to reconstruct the nuclear enhancement of gluon and sea quark parton 
distributions at moderate x. 

For the gluon and sea quark channels the corrections due to the ^-dependence of the 
valence quark distributions (which are anyway rather small for the x-range we discuss) 
enter only through Q 2 -evolution and are very small, see e.g., [25, 



2.1.3 Numerical predictions for nuclear PDFs and F^ 

In this subsection, we briefly discuss the DPDFs that we used in our analysis (we defer 
a more detailed specification to appendix [A]). We then present predictions for nuclear 
shadowing, both on the parton level (for both sea quarks and gluons) and for the nuclear 
structure functions, F- , for a range of Q 2 and A. 
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High precision data on diffractive processes, taken at HERA, and their QCD anal- 
ysis allow the diffractive parton distributions of the proton to be determined with good 
accuracy. In particular, both the ZEUS [29, 30, |3l| and HI [32| collaborations confirmed 



that a successful fit to the experimental data required a large gluon contribution (the 
same conclusion is drawn from an analysis based on the latest HI data [33] on inclusive 
diffraction). One should especially note that the recent HI data on diffractive photo- and 
electroproduction of dijets are dominated by the gluon diffractive structure functions 



We consider four different parameterizations of f® N '- fit D of Alvero, Collins, Terron and 
Whitmore [17]; fit D of [17], which has been improved to include a low-/? contribution; the 
parameterization from the theoretical light-cone model of Hautmann, Kunszt and Soper 



1 35], and the HI fit |Tj| to their own data [32]. Details of each of the parameterizations 
can be found in appendix [A]. 

Predictions for nuclear parton distributions fj/A made with these parameterizations 
of fj/N are fairly consistent with each other. The spread of the predictions reflects the 
theoretical uncertainty in determining fj/A 111 terms of f^j N due to experimental errors, 
certain differences between HI and ZEUS results, and the limited x and xjp ranges of the 
data. Furthermore, [|l7]] did not include the most recent ZEUS and HI data, while [35] 



actually did not perform a detailed fit to the data. In particular, the recent HI data p3| , 
which are most sensitive to the gluon diffractive structure functions, are best described 
by the HI fits [ll| to their earlier data, while the fit of Alvero et al. [17] somewhat 
overestimates the dijet production rate in the kinematics where the gluon DPDF at large 
(3 gives the dominant contribution 8 . 

For nuclear shadowing of the gluon distribution in nuclei, there is also an important 
uncertainty associated with the unknown i-slope of the diffractive gluon distribution of the 
nucleon, fg°/ N - A detailed analysis, presented in appendix |A|, demonstrates that Eqs. (|2.2| , 



2^4] , |2.7| ) require the diffractive parton distribution evaluated at t « 0. This implies that the 



answer for the shadowing correction is proportional to the t-slope of this distribution. The 
t-slope for the quark-dominated diffractive distributions in the proton, B q , was measured 
in inclusive diffraction by ZEUS |30| to be B q = 7.2 ± 1.1 GeV -2 (see appendix [A] for 
explanations on our use of B q ). No experimental information about the i-slope of the 
gluon diffractive distribution, B g , is available. In most of the phenomenological fits it is 
assumed that B g = B q . Since the amount of shadowing in the gluon channel increases with 
an increase of B g (for fixed total gluon-induced diffraction) we take a conservative attitude 
and consider a minimum value for the slope, which originates solely from the coupling of a 
small object to a nucleon via two gluons. To estimate this minimal slope we used the recent 
data on the diffractive electroproduction of J ftp by HI |3(| which can be parameterized in 
the following x-dependent form 



B g = (Bp + 2a' j ^ ln(10- 3 /^)) GeV" 



8 This is despite the fact that it has less shadowing. This effect is a peculiarity of the choices made by 
the various groups when choosing parameterizations for the input diffractive quark and gluon distributions 
in the regions of large (f3 ~ 1) and small /3 and reflects the residual uncertainties in the DPDFs. A larger 
DPDF at large f3 leads to a faster onset of nuclear shadowing at x ~ 10~ 2 , while the small-/? region 
determines the amount of shadowing at x < 1CF 4 . 
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with Bq = 4.5 GeV -2 and ol'ji^ = 0.125 GeV -2 . The slow x-dependence describes an 
increase of the slope with increasing energy ("shrinkage"). In order to reflect the uncer- 
tainty associated with the slope B g , we also considered a scenario of nuclear shadowing 
with B g = B q . For this choice, relative to that of Eq. ( |2.8| ), the shadowing of gluons, given 
by 1 — gA./(Ag^f), increases by a factor of 1.6 at x = 10~ 3 , and of 1.4 at x = 10 -5 , for light 
nuclei, and by a factor of 1.4 at x = 10 -3 , and of 1.2 at x = 10 -5 , for heavy nuclei. 

It is important to note that all of the parameterizations of the gluon DPDF in the 
region of x ~ 10 -4 — 10~ 5 reach the unitarity limit for Q 2 ~ 4 GeV 2 (see appendix |A|). 
This means that the probability of diffraction in gluon-induced DIS processes exceeds one 
half for these Q 2 . Hence, we are forced to tame the gluon DPDF at very small x in order to 
comply with the unitarity constraints. Note that, in the black-body limit, the cross section 
fluctuations are strongly suppressed: in principle reduces uncertainties in the predictions 
of shadowing for these x, within the framework of our calculation 9 . 

Figs. H] and || show our predictions for the ratio fj/A/i^-fj/N) for the up-quark sea 
and gluon distributions, respectively, in a range of nuclei. They clearly illustrate that the 
effect of nuclear shadowing for gluons is larger than for the quarks, especially at the input 
scale Q = 2 GeV. For instance, using the HI parameterization, one finds at Q = 2 GeV 
and x = 1CP 3 that the gluons are shadowed more than the quarks by a factor of 1.42 
for 12 C and by a factor of 1.24 for 208 Pb. At x = 10 -5 , due to the increase of B q with 
energy and the unitarity taming of the gluon mentioned previously, shadowing for gluons 
and quarks are very similar. When the ACWT parameterization is employed, the gluons 
are shadowed more than the quarks by a factor of 1.70 for 12 C and by a factor of 1.39 for 
208 Pb at x = 10" 3 . At x = 10" 5 , the factors are 1.29 for 12 C and 1.11 for 208 Pb. 

Of the four parameterizations of diffractive parton distributions discussed in appendix [A], 
we demonstrate only the two representative examples which give the largest (HI and 
smallest (ACWT [|l7|]) nuclear shadowing effects at x = 10 -5 . This illustrates the spread of 
theoretical predictions for fj/A/i^-fj/N) which arises from uncertainties in the diffractive 
parton distributions. The results are presented for Q = 2 GeV (solid curves), Q = 5 GeV 
(dashed curves) and Q = 10 GeV (dotted curves). At Q = 2 GeV, Eq. (2.7) was used. 
For the inclusive parton distributions of the nucleon, the recent CTEQ5M parameteriza- 
tion [37] was employed. We used the nuclear PDFs, fj/A-, at Q = 2 GeV as input to the 
DGLAP evolution equations in order to find fj/^ at Q = 5, 10 GeV. 

One can see from the figures that differences between the predictions are large for 
x ~ 0.01, where contribution of the large (3 > 0.3 — 0.6 region dominates. This (3 region 
was not strongly constrained by the older data from which the fits are derived. 

For the region of x ~ 10~ 3 predictions for sea quark shadowing exhibit minimal spread 
because the integrals over xjp in Eq. (2.7) are dominated by regions where the diffractive 
data is of high precision. The spread of predictions increases again for x < 10~ 4 reflecting 
the extrapolation of the DPDFs to the region of very high energies (beyond the HERA 
range). In the gluon case the spread is even larger, reflecting greater uncertainties in 
extraction of gluon diffractive distributions from the scaling violation of the F% data. 

9 However, since the approach to this limit is not understood and must involve new effects, our confidence 
in the predictions decreases. 
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Of two parameterizations presented here, only the HI parameterization fits the dijet data, 
which indicates that it may be more realistic. 
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Figure 4: The ratio of nuclear to nucleon sea up-quark distributions, ua/ (Aun), at Q = 2 GeV 
(solid line), Q = 5 GeV (dashed line), and Q — 10 GeV (dotted line). For each case, the scenarios 
with the largest (based on |l8|) and smallest (based on [jl7]]) nuclear shadowing for x = 1CP 4 are 
presented. 



It was observed in [10| that shadowing for gluons, calculated using the DPDFs of Alvero 
et al. [17], is much larger than for the sea quarks. The present analysis demonstrates that 
this is quite sensitive to the particular choice of the input diffractive parton distribution 
functions. However, it would appear from Figs. f| and |H| that for the smallest x values 
nuclear shadowing for sea quarks is not so much different from that for the gluons. There 
are several reasons for this. Firstly, and most importantly, our numerical calculations 
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Figure 5: The ratio of nuclear to nucleon gluon distributions, g^/tAgx) at Q = 2 GeV (solid line), 
Q = 5 GeV (dashed line), and Q = 10 GeV (dotted line). For each case, the scenarios with the 
largest (based on |Q) and smallest (based on |l7j]) nuclear shadowing for x — 10~ 4 are presented. 



were performed using the slope B g , given by Eq. (2.8), which is smaller than B g = B q , 
as implemented in [ 10 1 . Secondly, the effective cross section for the gluon channel, cfjp 



which defines the size of nuclear shadowing in Eqs. ( |2.2[ ) and Q2.4| ) for gluons, cannot be 
too large, otherwise the diffractive cross section becomes larger than allowed by unitarity 
of the scattering operator. Thus, in our model shadowing for the gluons reaches a plateau 
at Q ~ 2 GeV for x < 10 -4 , when the HI parameterization is used. In case of the ACWT 
fit, no taming of ofyy is necessary down to x = 10~ 5 . This shows that it is impossible 
to make realistic calculations of nuclear shadowing at small enough x, Q without invoking 
the unitarity restrictions. Also, the reverse statement is true: estimates of the kinematical 
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boundaries of the region in which unitarity restrictions play a role require that nuclear 
shadowing effects be properly taken into account! 

In order to demonstrate the importance of the value of B g , we also performed an 
analysis with larger B g : B g = B q . The results are presented in Fig. ^ for the lightest ( 12 C) 
and heaviest ( 208 Pb) nuclei, using the HI DPDF. The thin curves are those from Fig. || 
(using B g from Eq. ( |2.8|) ), while the thick curves are obtained using B g = B q . Again, 
the solid, dashed and dotted curves correspond to Q = 2,5,10 GeV. One observes that 
increasing the t-slope of the gluon DPDF increases the shadowing correction to gA, given 
by 1 — gA/ {A-9n)- For example, at x = 10 -3 and Q = 2,5, 10 GeV, 1 — gA/ {Ag^) increases 
by 60% for 12 C and 40% for 208 Pb. At x = 10~ 5 , the increase is about 40% for 12 C and 
20% for 208 Pb. 
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Figure 6: The ratio of nuclear to nucleon gluon distributions, gA /(Agw) at Q — 2 GeV (solid 
line), Q = 5 GeV (dashed line), and Q = 10 GeV (dotted line). The HI parameterizatio n |p~8[| for 
diffractive parton distributions is employed. The thin curves correspond to B g of Eq. (|2.8| ), the 
thick curves correspond to B g = B q . 

Using our results for the sea quark ratio qA/(AqN) and the gluon ratio gA/iAgw), one 
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can readily predict the ratio of the inclusive structure functions F 2 /(AF 2 ). In next to 
leading order accuracy, F 2 can be presented in the standard form pM] 

Ff{x, Q 2 ) = 4Wa{x, Q 2 ) + Q 2 )) (2-9) 



+ a s (Q 2 ) x f 1 dy 
2vr J x y 



Ci(x/y) ^ e 2 (q\(y, Q 2 ) + q A (y, Q 2 )) + (e 2 )C 9 2 (x/y)g A (y, Q 2 



Here is the electric charge of the quark with flavor i; (e 2 ) = e 2 )/nf, C| and Cf are 
the target-independent hard scattering coefficients 



2 3 



4rl + ^ 2 / 1 ,l-a\ 3\ 1 



ln( )-- +-(9 + 5x) 



1 — x V x 4/ 4 
C7| = n/ ((x 2 + (l-x) 2 )ln(i^)-l + 8x(l-x)) , (2.10) 

where nj is the number of active quark flavors and [. . .] + denotes the plus-regularization. 
In our analysis, rif = 4. In order to obtain F 2 , one needs to replace the nuclear parton 



distributions by the nucleon ones in Eq. (2.9). 



We considered two scenarios of nuclear shadowing, corresponding to ACWT and HI 
parameterizations of DPDFs. In the former case, since the DPDFs of the strange quarks 
is zero at the initial scale Qq = 2 GeV, only up and down quarks and gluons are shadowed 
at the initial scale. However, when the HI parameterization is used, f^ N = f^/N = fs/N 
at the initial scale so that up, down and strange quarks are shadowed. For both of the 
parameterizations of DPDFs we considered, the charm quarks are not affected by the 
nuclear medium at the initial scale, i.e., ca = Acm- However, as a result of the QCD 
evolution, gluons convert into ss and cc pairs, which leads to a significant deviation of the 
ratios sa/(Asn) and ca/(Acn) from unity (see Fig. |9|). 

The results of the calculation using Eq. ( |2.9D are presented in Fig. ^. In this analysis 
we assumed that valence quarks are shadowed the same as the sea quarks. A comparison of 
Figs. |I| and |7] reveals that while at Q = 5, 10 GeV, the ratios ua/ {Aun) and F 2 /(AF 2 N ) are 
quite similar, F 2 A /{AF 2 N ) is smaller than ua/(Aun) at Q = 2 GeV due to the important 
contribution of the unshadowed charm quarks. 

One can see from Fig. [?] that the solid curves corresponding to Q = 2 GeV run close 
to each other between x = 10 -4 and x = 10~ 2 . In order to understand this pattern, one 
should recall Eq. (^]) which demonstrates that 5F 2 is linearly proportional to F 2 ^ in the 
low nuclear density limit. Since both the ACWT and HI parametrizations give a good fit 
toif (4) ,it is no surprise that the resulting shadowing corrections to F 2 in the appropriate 
range of x are virtually the same for both fits. 

A peculiar feature of Fig. ^ is that the upper solid curves, which correspond to the 
calculation with the ACWT DPDF at Q = 2 GeV, in the small-x region lie above (i.e., cor- 
respond to smaller nuclear shadowing) all other curves corresponding to higher Q 2 . There 
are two reasons for this. Firstly and most importantly, the absence of nuclear shadowing 
for charm quarks at the initial scale significantly decreases the nuclear shadowing correc- 
tion to F 2 . Secondly, the next to leading order contribution to F 2 (terms proportional to 
a s (Q 2 ) in Eq. ( |2.9D ) has reduced shadowing as a result of the integration over y. 
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Figure 7: The ratio of nuclear to nucleon inclusive structure functions F% /(AF-f). The solid, 
dashed and dotted curves correspond to Q — 2,5,10 GeV and the scenarios of nuclear shadowing 
are as in Fig. ||. 



We would like to point out that there is no reason to assume that nuclear shadowing 
for the charm quarks is zero at the initial evolution scale. We propose the following simple 
model. At small values of Bjorken x charm quarks are mostly produced by QCD evolution 
because of the splitting g — * cc. Thus, nuclear shadowing for the charm quarks at some x 
and Q = Qo originates from nuclear shadowing of gluons at larger x and Q e ff- We choose 



c A (x,Q 2 ) _ 9A(2x,Q 2 e 
Ac N [x,Ql) Ag N (2x,Q 2 eff ) 



eff) (2.11) 



s effJ 

where Q^ff = 4m^ + Qq = 11 GeV 2 . Note that the gluon distributions are taken at twice 
the x-values of the charm quarks. 
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With this model for nuclear shadowing for charm quarks we re-evaluated the ratio 
F% /(AF2) using the HI parametrization of DPDFs. The results are presented in Fig. || 
as thin curves and compared to the corresponding curves from Fig. [7| (thick curves). One 
can see that by introducing nuclear shadowing for the charm quarks at the initial scale by 
Eq. ( 2.1 1] ) , one increases nuclear shadowing correction to F~ and makes it closer to the 
shadowing correction to ua- However, the size of nuclear shadowing introduced by our 
physically motivated model of Eq. ( 2.11] ) is not sufficient to have F^-f(AF^) « ua/{Aun) 
at small x. 
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Figure 8: The ratio of nuclear to nucleon inclusive structure functions F^/(AF^). The thick 
curves are those from Fig. ^, while the thin ones correspond to the calculation with shadowing for 
the charm quarks given by Eq. (2.11) at the initial scale Q = 2 GeV. The calculation is done using 
the HI parametrization of DPDFs. 



Both scenarios of nuclear shadowing for charm quarks at Q = 2 GeV and the results of 
their QCD evolution are presented in Fig. ^ using the HI fit to DPDFs. The thick curves 
represent our standard scenario without shadowing for the charm quarks at the initial scale, 
while the thin curves correspond to the calculation with shadowing for the charm quarks, 
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see Eq. ( 2.1 1| ) . The Q 2 -dependence of the curves is the same as in Fig. Two features 
of Fig. H are of interest: significant antishadowing and "wrong" Q 2 dependence of nuclear 
shadowing - nuclear shadowing increases as Q 2 increases! Both of these features can be 
understood as a consequence of the fact that at low x, charm quarks are predominantly 
produced by photon-gluon fusion j*g — > cc. 
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Figure 9: The ratio of nuclear to nucleon charm distribution functions, ca/{Acn)- The solid, 
dashed and dotted curves correspond to Q = 2,5,10 GeV and the scenarios of nuclear shadowing 
are as in Fig. 0. The calculation is done using the HI parametrization of DPDFs. 



The consistent inclusion of heavy quarks in input boundary conditions and in QCD 
evolution is known to be highly non-trivial, and involves a certain dependence on the scheme 
one is using for the number of active quark flavours (see |39| and references therein). This 
is one reason why such a low starting scale (i.e., below the charm mass) is often chosen, 
e.g., by CTEQ, so that heavy quarks are only generated explicitly via photon-gluon fusion 
at the input scale. Taken literally this implies that there can be no shadowing of charm 
at the input scale since there is no intrinsic charm at this scale ! We attempt to correct 
for this, in a phenomenological way, by using Eq. (|2Tll ) ab ove. However, a more accurate 
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treatment, along the lines of variable flavour number scheme, would be desirable but is 
beyond the scope of this paper. 



2.1.4 Predictions for central collisions 

The effect of nuclear shadowing becomes even more dramatic if one considers DIS from nu- 
clei at small impact parameters (i.e., "central" collisions). In this case, since the density of 
nucleons is larger in the centre of the nucleus than at its periphery, the number of scatterers 
effectively increases, which leads to an increase of nuclear shadowing. Experimentally, one 



can tag such small impact parameter events by measuring slow "knock-out" neutrons [40] 
which, in collider experiments with symmetric beam energies, can be detected with a very 
high efficiency (~ 100%) by low angle neutron calorimeters. 

Let us introduce a simple model for impact parameter-dependent nuclear PDFs, fj/ A (x, Q 2 -, b 2 )- 
By our definition, fj/ A (x,Q 2 ,b 2 ) are simply related to the usual, or impact parameter in- 
tegrated parton distributions, fj/A( x >Q 2 )i y i a 

J d 2 bf j/A (x, Q 2 , b 2 ) = f j/A (x, Q 2 ) . (2.12) 

In common with fj/ A (x, Q 2 ), fj/ A ( x ,Q 2 ,b 2 ) can be presented as a difference between the 
impulse approximation and shadowing contributions: 

f j/A (x, Q 2 ,b 2 ) =(j°° dz PA (b, z)\ Af j/N (x, Q 2 ) - Sf j/A (x, Q 2 , b 2 ) . (2.13) 

Here the factor dzp A (b, z) guarantees the correct normalization of fj/ A { x , Q 2 , b 2 ) (see 
Eq. ( 2.12Q ). The shadowing correction Sfj/ A (x, Q 2 , b 2 ) can be readily found from Eq. 



by removing the integration over b: 

Sf j/A ( x ,Q 2 ,b 2 ) = A{A ~ l) 16TrRe 



roc roo fip.o 

dz\ \ dZ2 \ dxjp X 



Z\ J X 



1 + T] 2 

ffM Q\ x ip, 0) PA(b, zi) PA (b, Z2 ) e ^ m ^-^e~ {A/2){1 -^<ff # dzpA(z) 



(2.14) 



We used Eq. (|2.14j) in exactly the same way that we used Eq. ( [2.7D . Firstly, the nu- 
clear shadowing correction to the gluon and sea quark parton impact parameter-dependent 



parton distributions was evaluated using Eq. ( 2.14 ) at the initial scale Qq = 2 GeV. Sec- 
ondly, antishadowing was modelled so that gluons in nuclei carry the same momentum 
per nucleon as in the free nucleon. Thirdly, the shadowing and antishadowing calculations 
defined the input to the QCD evolution to higher Q 2 scales. The results of such an analy- 
sis for the gluon distributions at the zero impact parameter are presented in terms of the 
ratio g A (x,Q 2 ,b = 0)/(AT(0)g N ) in Fig. (T(0) = / dzp A (b = 0,*)). The deviation of 
this ratio from unity is an effect of nuclear shadowing at small x and antishadowing for 
0.03 < x < 0.2. The results of Fig. |l(] should be compared with those of Fig. ||. By selecting 
DIS events with low impact parameters one can significantly increase nuclear shadowing. 
For instance, for Q = 2,5,10 GeV and x = 10~ 3 and x = 10~ 5 , the increase is about 
55 - 75% for 12 C and 20 - 40% for 208 Pb. 
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Figure 10: The ratio of the nuclear gluon distribution at the zero impact parameter to that of 
the nucleon, g A (x,Q 2 ,b = 0)/(AT(0)g N ) at Q = 2 GeV (solid line), Q = 5 GeV (dashed line), and 
Q = 10 GeV (dotted line). The ACWT and HI Q parameterizations for diffractive parton 
distributions are employed. 



2.2 The eikonal approximation approach 

Another popular and frequently used approach to nuclear shadowing is the eikonal approx- 
imation. The essence of this method in DIS on nuclei is the assumption that the forward 
amplitude for the virtual photon interaction with a nucleus can be written as the proba- 
bility of the transition 7* — > qq (the virtual photon light-cone wavefunction) convoluted 
with an exponential factor describing the interaction of the qq pair with the nucleons of the 
nucleus. The exponential factor is obtained by summing an infinite series of terms propor- 
tional to powers of the total gg-nucleon scattering cross section (i.e., by eikonalizing, hence 
the name of the approximation). 

An eikonalized form for high energy scattering amplitudes was first obtained by Cheng 



and Wu [41|, for several processes in QED. In particular, these authors demonstrated 



that for processes such as Delbruck scattering (photon-nucleus scattering in the static 
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approximation) the sum of Feynman graphs of a certain type can be cast in the form of a 
convolution of the e + e~-component of the virtual photon wavefunction with the eikonalized 
amplitude for the e + e~-nucleon scattering. 

Two aspects of the QED derivation of are important to mention. Firstly, it was ex- 
plicitly demonstrated that higher order diagrams, involving closed electron loops or graphs, 
in which a photon is allowed to be emitted and absorbed by the same electron, do not ex- 
ponentiate. Secondly, the transverse diameter of the e + e~-system emitted by the photon 
remains unchanged by the interaction. 

At ultra-relativistic energies, the eikonal approximation has been successful in the 
description of hadron-hadron scattering (for instance, the droplet model of |42|]) as well as 



hadron- nucleus scattering (see, for example, [43] and [44]). The approximation can be also 
applied to nucleus- nucleus scattering at high energies. 

The use of the eikonal approximation in high energy DIS on nuclei, within a dipole 
model for the wavefunction of a virtual photon, constitutes an extension of its previous use 
in hadron- nucleus scattering and QED. As explained above, the forward virtual photon- 
nucleus amplitude is given by the square of the photon light-cone wavefunction multiplied 
by an eikonal factor arising from the exponentiation of the total gg-nucleon scattering 
cross section. Such a form is based on the strong assumption that the gg-nucleon inter- 
action leaves certain characteristics of the qq-system, such as its transverse diameter and 
momentum fraction sharing variable, unchanged, i.e., the interaction is diagonal in the 
appropriate variables. In other words, within the eikonal approximation the projectile is 
said to be frozen in its prepared state for the lifetime of the interaction. 

A diagram typical of the eikonal approximation, as applied to nuclear shadowing, is 
shown in Fig. |H|. It represents the imaginary part of the forward virtual photon-nucleus 
scattering amplitude, in which the interaction with only a pair of nucleons is taken into 

A(2) 

account. Such a graph gives rise to the shadowing correction, 5F 2 , and should be 
compared to Fig. [I] of the leading twist approach. The differences between the two figures 
reveal conceptual differences between the approaches: while in Fig [l] the virtual photon 
interacts by dissociating into a multitude of diffractive intermediate states, the qq-dipole 
of the virtual photon in Fig. 11 interacts elastically through the exchange of an interacting 
gluon ladder (illustrated by the exchange of two gluons). 

In quantum mechanical processes, and high-energy hadron-nucleus scattering, the ap- 
proximation that the dominant fluctuations of the projectile can be considered to be frozen, 
can be justified by a suitable choice of basis states which diagonalize the scattering op- 
erator. However, this is not the case in QCD, where we are already given the basis, i.e., 
the quark and gluon degrees of freedom. In QCD, the gg-fluctuation of the virtual photon 
necessarily emits gluons, which can emit further partons, etc. The Fock components of 
the virtual photon, \qq), \qqg), • • •, intermingle before and during the interaction with the 
target, i.e., those states are not eigenstates of the scattering operator and cannot in general 
be considered to be frozen. Intuitively, the observation that emitting extra gluons would 
lead to the breakdown of the eikonal approximation is similar to the finding of [ 41 1 that the 
graphs where a photon is emitted and absorbed by the same fermion line, do not eikonalize. 

Since the eikonal approximation is so successful in describing cross sections of nucleon- 



- 23 - 




Figure 11: Feynman diagram giving rise to nuclear shadowing within the eikonal approximation. 



nucleus scattering, should it not also be a good method to describe lepton-nucleus DIS? The 
above discussion shows that the eikonal approximation is justified only for processes where 
the presence of the | (/(/(^-component of the virtual photon is unimportant. An example of 
a relevant observable is F 2 A at Q 2 of only a few GeV 2 . On the other hand, for reactions 
sensitive to the \qqg) -component, the eikonal approximation is expected to fail and give 
wrong size and Q 2 -dependence of nuclear shadowing. One example of an observable for 
which this effect is dramatic is the longitudinal structure function F A . 

Within the eikonal approximation, the shadowing correction to the nuclear inclusive 
structure function F 2 A , 5F A , can be written in the form 



6F 2 A (x,Q 2 



Q- 



4vr 2 



-Re 



"cm 

DO 



J da d 2 d t Y^ l 1 ^"'*? 2 '^' 



2 ^ 



d 2 b dz x dz 2 (l-ir,) 2 aSv(x,4,m 4 ) p A (b, Zl ) p A (b, z 2 )e l2xm ^ z ^ x 



-(A/2)(l-ir,y q f N (x,dl,mi) dz PA {z) 



(2.15) 



Here a em is the fine-structure constant; a is the fraction of the photon's longitudinal 
momentum carried by q or q; dt is the transverse diameter of the (/(/-system; mi is the mass 
of the quark with flavor i; p A is the nuclear one-body density. The square of the light-cone 
wavefunctions of the virtual photon is given by the standard expression 



Ma,Q 2 ,dl m )\ 2 ~ 6ae 



7T- 1 



Q 2 a 2 {l-af + 1 ^f)K 2 {e l dt) 



+ - [a z + {l-aY) 4 K((eid t) 



(2.16) 



where Kq and K\ are the modified Hankel functions; e« = Q a(l — a) + mf; t] is the ratio 
of the real to imaginary parts of gg-nucleon scattering amplitude. In our analysis, we set 
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rj = 0.25, which is consistent both with our leading twist calculations and vector meson 
dominance ideas, see e.g., [15|. Note that since F2 measures the virtual photon-target cross 
section averaged over helicities of the incoming photon, averaging over the photon helicities 
is assumed in Eqs. (2.15) and ( 2.1 6| ) . Following our analysis in [45], we include four flavors 
of quarks and assume m u = = m s = 300 MeV and m c = 1.5 GeV. 

The use of the optical theorem for the elementary qq-nucleon amplitudes enables one 
to express the answer in Eq. ( p. 15 ) through the total qq-nucleoii cross section, c^v- This 
cross section plays a key role in the so-called dipole formalism of DIS, and we refer the 
reader to [46] for a brief review of various formulations of the dipole formalism existing 
in the literature. In this analysis, we used a QCD-motivated model for o l °^ N developed in 
1 45]. For small dipoles, the total cross section is predominantly inelastic and, as such, is 
governed by the gluon distribution in the proton. In the non-perturbative region of large 
dipole sizes, we model by requiring its equivalence to the soft pion-nucleon total cross 
section. 



In exactly the same or very similar form, Eq. ( |2.15| ) was used in a number of papers 



EH, 



I ^9|, [5(], |52]] . In the derivation of Eq. ( 2.15| ) one assumes that the invariant 
mass of all gg-dipole intermediate states, which is inversely proportional to the dipole 
diameter dt, is the same and approximately equal to Q (i.e., the diffractive variable = 
Q 2 /{Q 2 + M 2 ) ps 0.5, which implies x P ra 2s). This means that the non-zero longitudinal 
momentum transfer to each individual nucleon of the nuclear target is of the order of 2xniN, 
which explains the factor exp(i2xmAr(zi — z?)) in Eq. ( |2.15| ). 

We would like to stress the following subtle point, which should be kept in mind while 
applying Eq. (2.15) and which was ig nored in g|, ||, ||, ||, ||]. The total cross section 
(rt AT receives contributions from both elastic, crfi^r, and inelastic, a"^ e A \, cross sections. For 



qqN 



qqN ' 



'qqN' 



small dipole sizes and not very small x, cr l ^ N is relatively small and can be calculated using 
methods of perturbative QCD. Moreover, in this case o e ^ N is negligibly small and thus one 



'qqN With OtfN- 



can identify <r*|* 

large, it is incorrect to ignore the contribution of a el 



In circumstances when a l ^ N and cr^ N become sufficiently 



qqN 



to the total cross section. For 



example, and <^^ N become equal in the limit when the qq-micleon scattering amplitude 
reaches its limiting value allowed by unitarity of the scattering operator. 



In order to illustrate how the eikonal approximation works, we apply Eq. ( 2.15 ) in 
the most straightforward way by assuming it is valid at all Q 2 . Of course, as explained 
previously, this becomes progressively less justified as Q 2 increases because of the mixing 
of \qq), \qqg) and \qqgg ■ ■ ■) fluctuations of the virtual photon (in practice this restricts 
legitimate values of Q 2 to a few GeV 2 ). 

As it stands, Eq. (2.15) overestimates nuclear shadowing at the higher end of the 
shadowing region, 0.01 < x < 0.07. In particular, 5F^ does not vanish at x = 0.1, as 
happens in the leading twist approximation. Thus, in order to obtain sensible results with 
Eq. ( 2.15j ), it is applied only for 10~ 5 < x < 0.01. For the interval 0.01 < x < 0.1, we 
assumed that SF^ decreases linearly and becomes zero at x = 0.1. 

The ratio F^/(AFf), calculated using Eq. (|2.15| ), is presented in Fig. ^ for 12 C and 
208 Pb. The Q 2 -dependence of this ratio is shown by solid (Q = 2 GeV), dashed (Q = 5 
GeV) and dotted (Q = 10 GeV) curves. A comparison of Figs. |l2| and ^ reveals the 
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The ratio of nuclear to nucleon inclusive structure functions F 2 



A /(AFf) calculated 



Figure 12: 

within the eikonal approximation using Eq. (2.15). The Q 2 -dependence is given by the solid (Q = 2 
GeV), dashed (Q = 5 GeV) and dotted (Q = 10 GeV) curves. 



following characteristic trends (from the considered scenarios of nuclear shadowing within 
the leading twist approach, we choose the one corresponding to the HI parameterization). 
At Q = 2 GeV, and for the 12 C and 208 Pb nuclei, both leading twist and eikonal approaches 
give similar (within about 20%) predictions for F^/^AF^). As Q 2 increases, for light and 
heavy nuclei the leading twist approximation predicts much larger shadowing than the 
eikonal approach. In particular, at Q = 10 GeV and x = 10~ 3 , F% /(AF^) is shadowed 
more in the leading twist approach than in the eikonal approximation by 74% for 12 C 
and by 71% for 208 Pb. At x = 10 -5 , the corresponding increase of nuclear shadowing is 
similar for 12 C and more modest for 208 Pb: it is 72% for 12 C and 49% for 208 Pb. The 
Q 2 -behaviour follows the expected pattern: while the Q 2 -behaviour of F^ /(AF^) within 
the leading twist approach is governed by the QCD evolution equation, and is therefore 
logarithmic, within the eikonal approximation it decreases with increasing Q 2 much faster 
and is dictated largely by the Q 2 -dependence of the virtual photon light-cone wavefunction. 
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2.3 Differences between the leading twist and eikonal approaches 

In this subsection, we summarize what has been said so far about the differences between 
the leading twist and eikonal approaches to nuclear shadowing in DIS on nuclei as well as 
discuss problems with an unambiguous implementation of the eikonal approximation. 

The two key differences, which make the approaches so distinct, are related. They are: 

• a different space-time evolution of the scattering process; 

• the neglect of the \qqg) component (and higher Fock states) of the virtual photon in 
the eikonal approximation. 

The eikonal approach is reliable only within the framework of non-relativistic quantum 
mechanics, where the number of interacting particles is conserved during collisions. In this 
case, the approximation that the Fock states of the incoming high-energy virtual photon 
can be considered frozen, is justified, and the procedure of eikonalization can be successfully 
carried out. 

In contrast, in a quantum field theory such as QCD, the number of bare particles is 
not conserved. In other words, the number of effective degrees of freedom, or relevant 
Fock states, in the photon wavefunction depends on x and Q 2 . For example, the interact- 
ing \qq) Fock state radiates gluons, thus creating and mixing with \qqg), \qqg ■ ■ ■ g) states. 
This mixing is properly taken into account by the QCD evolution in the leading twist ap- 
proximation. Of course QED is also a quantum field theory, however, the electromagnetic 
coupling is much weaker than the strong coupling of QCD, and there is no self- interaction 
of photons, so the corrections which spoil the eikonal approximation for QCD are corre- 
spondingly much smaller in the QED case. 

One immediate consequence of the contrasting space-time evolution pictures within the 
leading twist and eikonal approaches is the Q 2 -dependence of nuclear shadowing. As Q 2 
increases, the Fock components of the virtual photon with an increasing number of gluons, 
\qqg---g), become important for nuclear shadowing. This follows straightforwardly from 
the connection between nuclear shadowing and gluon-dominated diffraction, as found in 
ZEUS and HI experiments at HERA. As a result, using the factorization theorem, the 
Q 2 -dependence of nuclear shadowing is governed by the DGLAP evolution equation within 
the leading twist approach. 

In the variant of the eikonal approximation that we considered, the \qqg) Fock state 
of the virtual photon is absent. Hence, even if one adjusts the |qg)-nucleon cross section 
in order to reproduce correctly the nuclear structure function at the initial scale Qo, 
the approach would necessarily underestimate nuclear shadowing at large Q 2 since the 
eikonal approximation is not based on QCD evolution. This was already demonstrated by 
comparing Figs. 12 and [?]. As we will illustrate in the next section, this effect is even more 



dramatic for the longitudinal structure function F£. 

The naive formulation of the eikonal model, which only contains the gg-component of 
the virtual photon wavefunction, underestimates the amount of nuclear shadowing since it 
neglects diffractively produced inelastic states, such as qqg, qqgg, etc. The inclusion of a 
qqg component, as was done in the analysis of diffraction on the nucleon in |j53|l, is only 
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a step in the right direction. In order to reproduce the correct Q 2 -behaviour of nuclear 
shadowing, which is governed by the DGLAP evolution equation, one should include the 
complete set of Fock states, i.e., an infinite series of components including infinitely many 
constituents. 

There are several subtle points and technical problems with implementation of the 
eikonal approximation. Firstly, one has to be careful to use the total gg-nucleon cross 



section in Eq. (2.15). In the kinematics, where the elastic and inelastic qq-nucleon cross 
sections are compatible, the use of alone would significantly underestimate nuclear 
shadowing. 

Secondly, in order to reproduce correctly nuclear shadowing at the higher end of shad- 
owing region, 0.01 < x < 0.07, one needs to take into account the non-zero longitudinal 
momentum transfer to the nucleus through the factor exp(i2xm.N(zi — 22)). In order to 
arrive at this factor in the eikonal approximation, one needs to make a strong assumption 
that all essential Fock states of the virtual photon have the same invariant mass of the 
order of Q. 



Thirdly and very importantly, there is no unambiguous way to generalize Eq. ( 2.15 ) for 
nuclear structure functions to something similar to Eq. ( |2.4D for individual nuclear parton 
distributions. Of course, one can attempt to replace the ratio F^-/(AFJy) in Eq. ( 2.15j ) by 



fj/A/ (AfjMr) but what can one use for cross section (T to ^ N ? The eikonal approximation 



gives no clue as to cr t( £ N for different flavors of partons. Moreover, since such a picture would 



I j /Al \ A J j /N ) UUL wiiai can une use rur crusts secuoii o ^ 
T tot 

T qqN 

not be based on the factorization theorem and QCD evolution, the scaling violations of 
quark and gluon nuclear distributions would not be consistent with each other and the 
DGLAP equations. 

To conclude we would like to give a clear answer to the question as to why the eikonal 
approximation works so well for hadron-nucleus processes but fails for lepton-nucleus DIS 
processes. For high energies the hadronic projectile fluctuates into configurations which 
subsequently interact with the target with similar cross sections. This means that the 
distribution over cross section fluctuations is rather narrow. As a result, in hadron-nucleon 
scattering inelastic diffraction is a small correction to elastic scattering. When these ideas 
are applied to hadron-nucleus scattering, one sees that inelastic intermediate states (corre- 
sponding to inelastic diffraction) give a small contribution compared to elastic intermediate 
states (corresponding to elastic scattering). Hence, the eikonal approximation works well. 

In lepton-nucleus DIS the situation is very different. Firstly, both small- and large-size 
fluctuations of the virtual photon contribute to DIS on the nucleoli, which means that 
cross section fluctuations are very significant. Secondly, for small-size fluctuations, the 
situation for DIS off a nucleon is opposite to that assumed in the eikonal approximation: 
in the former the cross section for inelastic diffraction (a leading twist observable) exceeds 10 



The reason for the dominance of inelastic diffraction with an increase of Q at fixed small x can be 
traced to the definition of elastic cross section, which was introduced using t-channel factorization so that the 
virtual photon-target cross section has the form of the convolution of the photon (qq) wavefunction with the 
qg-target cross section. However, as gluons are attached to the qq pair at large Q, they are predominantly 
emitted well before the photon wave packet has approached the target. These configurations containing 
gluons, which on average have a large transverse size and may be in the color octet dipole state rather than 



- 28 - 



the cross section for elastic scattering (a sub-leading twist observable). For the latter the 
reverse is true. 

As a result, inelastic intermediate states dominate in lepton- nucleus DIS initiated by 
small dipoles. Of course, the presence of large-size dipoles in the virtual photon wavefunc- 
tion introduces a certain degree of similarity between hadron-nucleus and lepton-nucleus 
scattering, but this is not sufficient to justify the application of the eikonal approximation 
to DIS on nuclei. 

3. Nuclear shadowing of the longitudinal structure function Ff 

As discussed above, the crucial difference between the leading twist and eikonal approaches 
is contrasting space-time evolution pictures of the scattering process. As a consequence, the 
higher Fock components of the virtual photon, containing gluons, are effectively included 
in the leading twist approach and are neglected in the eikonal approximation. Therefore 
any observable which is sensitive to the gluon distribution in the nucleus should be a 
good tool to distinguish between the leading twist and eikonal approximations. We have 
already demonstrated this using the inclusive structure function, F^, at large Q 2 . An even 
more striking example is given by the longitudinal nuclear structure function, F£, which 
is measured by DIS of longitudinally-polarized photons on nuclei. Other relevant processes 
include exclusive electroproduction of p (dominated by longitudinally-polarized photons) 
and J ftp mesons off nuclei. 

The nuclear structure function can be obtained by a straightforward generalization 
of the one-loop perturbative QCD result for the nucleon pi] : 

f£(x,Q*) = 2 -^± f d A-?(^4{l- X -)y gA {y,Q 2 ) + \Fi{y,Q 2 )) . (3.1) 

71 Jx y y v ^ y 3 J 

Here nj is the number of active quark flavors. Replacing qa by gpj and F% by F^ , one can 
present the ratio of the nuclear to nucleon longitudinal structure functions, F- / '(AF^f), in 
the form 

F^Q 2 ) _ £^(f) 2 (Er4e?(l-f)^A(i/,Q 2 ) + li^(!/,Q 2 )) 
AF^(x, Q 2 ) £ ^£(£)2 e?(l - %)yg N (y, Q 2 ) + \F^{y, Q 2 )) ' 

The ratio f£/(AF£), calculated using Eq. (|3.2[) , is presented in Fig. |l3| (curves are 
labelled as in Fig. |||). By comparing Fig. [L3] to Fig. |B| one sees that the size and Q 2 - 
dependence of nuclear shadowing in the ratio F^ /(AF^) are similar to those in the ratio 
gA/(Agiy). Slightly smaller shadowing for F^/(AF^) is an effect of the convolution and 
sea quarks present in Eq. (|3^). 

in the color triplet dipole state, much more readily rescatter diffractively (i.e., with a small momentum 
transfer to the nucleon). As a result, the genuine total cross section of the interaction, as measured through 
double rescattering, turns out to be much larger at large Q 2 than the cross section defined via the elastic 
qqN cross section. 
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Figure 13: The ratio of nuclear to nucleoli longitudinal structure functions F% /(AF^) calculated 



within the leading twist approximation using Eq. (3.2). Different curves correspond to different 
values of Q 2 and scenarios of shadowing as in Fig. |. 



Within the eikonal approximation, one cannot use Eq. ( |3.2| ) because the eikonal ap- 
proximation does not present an unambiguous scheme to evaluate g A - Instead of Eq. (|3.2|), 
the longitudinal structure function in the eikonal approximation can be obtained using 
Eq. (|2l5|) by repl acing the helicity-averaged photon wavefunction \& by the longitudinally- 
polarized one ^ l'- 

6a P 



\^ L (a,Q z ,d 2 ,m 



Y,^Q 2 a 2 (l-a) 2 Kl(e t d t ) . 



(3.3) 



Then the ratio F£/(AF^) can be presented in the form 
F£(x,Q 2 ) .1-1 
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2 



-Re 



dad 2 d t ^2\^ L (a,Q 2 ,d 2 t )\ 2 J d 2 b J 



oo poo 

dz\ \ dz% x 

OO J Z\ 



- 30 - 



( J dad 2 d t ^2\^ L (a,Q 2 ,d 2 ,m i )\ 2 a t q f N (x,d 2 L ,m i ) 



(3.4) 



Results for the ratio of the longitudinal structure functions in the eikonal approxima- 
tion are given in Fig. 14. As for FA/^AFN) within the eikonal approximation, we assumed 
that F^ /(AF^) increases linearly for x > 0.01 and becomes unity at x = 0.1. 




Figure 14: The ratio of nuclear to nucleon longitudinal structure functions F 2 L /(AF^) calculated 
within the eikonal approximation using Eq. (^J). The Q 2 -dependence is shown by the solid (Q=2 
GeV), dashed (Q=5 GeV) and dotted (Q=W GeV) curves. 



The comparison of Figs. ^ and 14 clearly reveals dramatic differences between the 
leading twist and eikonal approach predictions. Firstly, the leading twist approach already 
predicts much larger shadowing at the initial scale Qq = 2 GeV. Secondly, nuclear shadow- 
ing dies out with increasing Q 2 much faster within the eikonal approximation. In order to 
quantify this discussion, our results for the nuclear shadowing correction to the longitudinal 
structure function, 1 — F^/(AF^), of 12 C and 208 Pb are presented in Table 1 at x = 10~ 3 
and x = 10 -5 . For the column containing the leading twist results, first values correspond 
to the ACWT parameterization (and the values in parenthesis to the HI parameterization). 
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Table 1: Nuclear shadowing correction to the longitudinal structure function within the leading 
twist and eikonal approaches. 



One can see from Table [l] that, depending on x and A, the shadowing correction to in 
the leading twist approach is larger, compared to the eikonal approximation, by a factor of 
1.3 — 2 at Q = 2 GeV and by a factor of two to seven at Q = 10 GeV. Hence, measurements 
of F£ should be an excellent means to distinguish between the leading twist and eikonal 
approaches to nuclear shadowing. 

4. Conclusions and discussion 

We compare two frequently used approaches to nuclear shadowing in DIS from nuclei: the 
leading twist and eikonal approaches. Our comparison is based on the observation that one 
of the foundation principles of the two approaches, the space-time picture of the interac- 
tion, is quite different. The leading twist approach is based on the observation that at high 
energies the projectile interacts simultaneously with several nucleons of the nuclear target. 
This leads to the reference frame independent connection between nuclear shadowing in 
DIS from nuclei and DIS diffraction from the nucleon. The QCD factorization theorems for 
inclusive and diffractive scattering enable one to express the shadowing correction to the 
individual parton distribution functions in nuclei in terms of the corresponding diffractive 
parton distribution functions of the proton. Moreover, another consequence of the factor- 
ization theorems is that the scaling violations of nuclear parton distribution functions are 
given by the leading twist DGLAP equation. 

Unlike the leading twist approach, the eikonal approximation is a frame-dependent 
method, which applies in the target rest frame. In this reference frame, the incoming virtual 
photon first fluctuates into partonic components (\qq), \qqg), ■ ■ ■), which then interact with 
the nucleus. We considered the often- used approximation in which only the \qq) fluctuation 
is considered. In this case, the qq effective dipole interacts successively and elastically with 
nucleons of the target. Apart from this, an additional assumption is made that such 
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interactions can be presented in an exponential or eikonalized form. It is important to 
realise that the eikonal approximation implies that only elastic (\qq)) intermediate states 
contribute to the virtual photon-nucleus cross section. In addition, the method is not 
based on the QCD factorization theorem and therefore does not include the proper QCD 
evolution. 

In this work, we performed several numerical studies. 

Firstly, within the leading twist approximation, the ratio of the nuclear to nucleon 
gluon and sea quark distribution functions was examined as a function of Bjorken x for 
several values of Q and for a wide range of nuclei. By testing four distinct parameterizations 
of diffractive partons distributions of the proton, we extended the earlier analysis of [10] 
and confirmed that the gluons in nuclei are shadowed more than the sea quarks. We found 
that the difference between gluon and sea quark nuclear shadowing was not quite as large 
as that found in juj (compare Fig. ^ to Fig. 5 of |l(J which uses the ACWT fit). There are 
two major reasons for this. Firstly the Hl-type fit corresponds to a smaller diffractive gluon 
density at large (3 than ACWT. Secondly, in the current analysis we allowed for a possible 
difference of the slopes for the diffractive gluon and quark distributions (see appendix [A| for 
details). In addition, we made predictions for shadowing of the gluon parton distributions 
at central impact parameters as well as for the charm parton distributions. 

Secondly, we compared predictions for the nuclear structure functions, F^, within the 
leading twist and eikonal approximations. One should emphasize two aspects of this com- 
parison: the size and Q 2 -dependence of the nuclear shadowing correction. At the initial 
scale Q = 2 GeV, both leading twist and eikonal approaches give similar predictions. Nev- 
ertheless, the Q 2 -dependence of F^ gives good discriminating power between the leading 
twist and eikonal approaches which reflects their conceptual differences: a logarithmic de- 
crease of nuclear shadowing in the leading twist approach is much slower than that predicted 
using the eikonal approximation. For instance, at Q = 10 GeV and x = 10 -3 , the leading 
twist approach predicts 71% more shadowing for F^ 0Sph than the eikonal approximation. 

Thirdly, given the fact that the eikonal approximation omits contributions of the vir- 
tual photon fluctuations containing gluons (\qq~g), ■ ■ ■) it must predict very small nuclear 
shadowing for observables sensitive to the gluon distributions of nuclei. Examples of such 
observables include cross sections for longitudinally polarized p meson electroproduction, 
J/ip electroproduction, and the longitudinal structure function F^. We considered the 
latter in detail within the leading twist and eikonal approximations. The differences be- 
tween the predictions for the size and Q 2 -dependence of the shadowing correction to F£ 
are very dramatic: the eikonal approximation undershoots the leading twist prediction for 
nuclear shadowing by at least a factor of two for all of the nuclei and values of Q 2 which we 
considered. Therefore, since F£ can be measured at the future Electron-Ion Collider (3|, 
this observable should be able to discriminate unambiguously between the leading twist 
and eikonal approaches to nuclear shadowing in DIS from nuclei. 

Nuclear shadowing is closely related to another interesting small- x phenomenon namely 
the violation of the DGLAP approximation at sufficiently small values of x. This phe- 
nomenon is referred to by various names in the literature: parton saturation, parton tam- 
ing, unitarity constraints, etc. Let us consider DIS from a heavy nucleus at very large 



- 33 - 



energies, i.e., at very small x. For almost all essential impact parameters, the nuclear 
structure function is predicted to increase mildly with decreasing x, oc ln(l/a;), which 
is significantly slower than the behaviour predicted by the DGLAP equation (see e.g., p5(|). 
It is now understood that parton saturation and related phenomena occur in heavy nuclei 
at larger x than in the proton, which justifies the use of nuclear beams in attempt to study 
new aspects of small- x physics. On the other hand, the leading twist nuclear shadowing 
significantly reduces the density of partons (especially gluons) at small x and thus com- 
petes with the all- twist phenomenon of parton saturation. As a result, sensible studies of 
QCD at high parton densities can be carried out only if nuclear shadowing of partons is 
properly taken into account. 
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A. Parameterizations of diffractive parton distributions and a e ff 

Within the leading twist approach, the nuclear shadowing correction, Sfj/^, to nuclear 



PDFs is given by Eq. (2.7). It involves the effective cross section, o^**, which gives the 



e// ' 

strength of the interaction with any two nucleons of the nucleus leading to nuclear shadow- 
ing of the parton of flavor j. This may be expressed through the nucleon diffractive parton 
distribution function (DPDF) for a parton of the same flavour as follows [|l^, |2!| : 



^'^'' ^pWUHfl l ^W^CftflW-O), (A.1) 

f D(4) 



where fj/ N is the inclusive parton distribution of the nucleon; f^ N is the diffractive parton 
distribution of the nucleon (which is, strictly speaking, a proton); r\ is the ratio of the real 
to imaginary parts of the diffractive scattering amplitudes. In our analysis, rj = 0.22 for 
the amplitude corresponding to the ACWT and model 3 parameterizations and rj = 0.32 
when the HI parameterization is used. 

As explained in Sect. B the upper limit of integration, Xjp^, is different for the gluons 
(xiPfi = 0.03: to allow antishadowing) and for the sea quarks (xjp^ = 0.1). However, for 
the sea quarks, because of the factor exp(ix]pmN(zi — Z2)) in Eq. (|2.7D and the absence of 
antishadowing, the exact value of xjp^ is unimportant for our numerical analysis. 

Ignoring the small minimum momentum transfer, t m i n , always present in inelastic 
diffraction, the condition t ~ —k 2 = in Eq. ( |A.1| ) means that f^fpj is to be evaluated 
at t = 0. However, all the parameterizations to the diffractive data which we used are fits 
to i-integrated diffractive structure function F% • This means that we should assume 
a certain t-dependence of /f/jy '• For practical reasons, we assumed a simple exponential 
dependence 

f^\p,Q 2 , XlP ,t) = e- B ^f^\p,Q 2 , XlP ,t = 0) , (A.2) 
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where the slope B could depend on xp (for the sea quarks) as well as on Bjorken x (for 
the gluons). 

Introducing the diffractive distribution function /jyjy as 

f D ^(l3,Q 2 ,xjp) = f dtf^(P,Q 2 , X]P ,t) , (A.3) 

J — oo 

one readily obtains that 

f^ ) ((3,Q 2 ,x ]P ,t = 0) = BM/J(AQ 2 ,x P ) . (A.4) 

As a next step, we use the Regge factorization hypothesis which assumes that f is 
a product of a Pomeron flux factor, /jp(xjp), and the parton distribution functions of the 
Pomeron, fj/ip{(3,Q 2 ): 

f^\f3,Q 2 ,xjp) = Mxip)fj/ip(f3,Q 2 ) ■ (A.5) 



The comparison of Eqs. (A.4) and (A.5) gives 



ff /N (J3, Q z ,xjp, t = 0) = B(xjp) fjp{xjp) f j/F (/3, Q 2 ) . (A.6) 
Substituting Eq. ( |A.6| ) into Eq. ( |A1| ), we obtain our master equation for <7g// : 

^ ff(x ' Q2) = l x dxpB{x]p) Mxf) Q2) ■ (A - 7) 

The theoretical analysis of the HERA diffractive data, in terms of quark and gluon de- 
grees of freedom, with the additional assumption of Regge factorization and Pomeron domi- 
nance, effectively concerns itself with the parton distributions in the Pomeron, fj/jp((3, Q 2 )- 
One way to study fj/jp phenomenologically is to fit the experimental data by performing 
QCD evolution using a reasonable trial shape for fj/jp at the initial scale Qo- Models 1, 
2 and 4, discussed below, are examples of such a determination of fj/jp- Model 3 is a 
theoretical prediction for fj/jp, which also successfully describes certain diffractive data. 

It is important to emphasize that the four considered parameterizations of the gluon 
distribution function in the Pomeron, f g /]p, become so large at small values of Bjorken 
x (x < 3 x 10 ) that the gluon-induced diffractive DIS cross section exceeds half the 
gluon-induced inclusive DIS cross section. Clearly, this is prohibited by the unitarity of 
the S matrix (when the unitarity limit is reached, the diffractive and elastic cross sections 
are equal each other and are equal to half the total cross section). In terms of f g /jp, the 
unitarity restriction can be written in the form 



x p,o 1 

dxpf^x^fg/jpix/xp^Q 2 ) < -g(x,Q 2 ) , (Ai 



where g(x, Q 2 ) is the number density of gluons in the nucleon. So, whenever / g /jp becomes 
so large as to violate the unitarity limit of Eq. ( A.8| ), f g /jp is replaced by its limiting value 
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given by Eq. (A/7). This allows a 9 e ^ to continue to grow as x decreases but only at a 
reduced rate dictated by the growth of the inclusive gluon PDF. 

Predictions for a c g for the up-quark sea and gluon, obtained using Eq. ( |A.7| ) (with 
the unitarity restrictions on are presented in Figs. ^ and [l(J. All curves correspond 
to Q = 2 GeV. For the inclusive parton distributions, the CTEQ5M parameterization [37] 
was used. For the gluons, we tested two choices of the diffractive slope B g : the solid curves 
correspond to B g parametrised by Eq. (p.8|), while the dotted curves are obtained with B, 



of Eq. ( |A.12[) below. 



We now describe the parameterizations of DPDFs, which we used, in detail. 
Model 1 

Numerical predictions in are based on the parameterization by Alvero, Collins, 
Terron and Whitmore |jT7| . The choice of fit D gives the best fit to the data on diffractive 
DIS and diffractive photoproduction of jets taken at HERA by the ZEUS and HI collabo- 
rations. Note that this parameterization somewhat overestimates the most recent HI data 



on jet production [34] 
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Figure 16: The effective cross section crf ff for the gluons obtained using Eq. ( A. 7 ) (with the 
unitarity restriction) at Q = 2 GeV. The solid curves correspond to the slope B g given by Eq. ( |2.8| ) , 
while the dashed curves are calculated using B g of Eq. (A. 12). 



The following shapes of the quark (q = u = u = d = d) and gluon parton distribution 
functions of the Pomeron were obtained at the initial scale Qo = 2 GeV |T|] 

Pf q /M>(P, Ql) = 0.292 -13)- 0.159(1 - /3) 2 ) , 

Pfg/riP, Ql) = 9-7/9(1- /?) . (A.9) 

The strange and charm quark distributions are taken to be zero at the initial scale Qo an d 
are generated by QCD evolution for Q 2 > Qq- 

We would like to point out two important features of the parameterization given by 
Eq. ( |A.9p . Firstly, a successful DGLAP fit to the diffractive data requires that the gluon 
distribution is much larger than the quark distribution, in order to reproduce the observed 
scaling behaviour in Q 2 . This automatically implies that cf^j is much larger than f"// 
(e.g., compare Fig. 15 to Fig. 16 in the region around x = 10~ 3 ). Secondly, the quark 



distribution is rather hard and it even becomes negative at low values of (3. However, this 
does not lead to a paradox since the diffractive data with low (3 were not used in the fitting 
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procedure of jT 
constrained. 

The Pomeron flux used in Eq. (| 
form 



Thus, the low-/3 behaviour of the parameterization ( A.9 ) is not well 
) was parameterized in the Donnachie-Landshoff 



dt 



Ami - 2.8t 



4m 2 



1 



t v l-t/V 



l-2ajp(t) 



(A.10) 



where m p is the proton mass; /3q = 1.8 GeV -1 is the Pomeron-quark coupling; t\ = 0.7 
GeV 2 ; ajp{t) = ajp + olt is the Pomeron trajectory. The analysis of [17] showed that a 
successful fit the diffractive data favors ajp = 1.14 and a' = 0.25 GeV -2 . 

The slope B q (xp) for the quarks, which enters Eq. QA.7 ), was found from the t- 
dependence of the Pomeron flux (see Eq. ( A. 10) )) 



Bq{xjp) 



l-2ajp(t=0) 
J JP 



dt 



r4m 2 



2.8* 



1 



Am 2 



t M - t/ti 



l-2ajp(t) 
iP 



(A.ll) 



Thus defined B q {xp) is a slow function of xp: B q increases with decreasing xp. For a 
wide range of xjp, the value of B q is close to the slope of the diffractive structure function 
F 2 D(4) , B = 7.2 ±1.1 GeV" 2 , which was measured by the ZEUS collaboration @. 

For the gluon diffractive distribution, the slope B g is expected to be smaller than B q . 
The rationale for the use of a smaller slope for gluons in this and other models comes from 
two different lines of argument |5(| and [57|. In [^] it was argued that B g ~ 5 GeV -2 
leads to effective cross sections at Qo ~ 2 GeV similar to that which would be obtained by 
rescaling the quark-antiquark dipole cross sections we found in by the Casimir operator 
factor of 9/4. Also we pointed out that B g ~ 7 GeV~ 2 would lead to a e ff for gluons as 
large as 60 mb, which would imply blackness for the interaction up to very large impact 
parameters. In ]5^] , based on the calculation of the diffractive parton densities, it was 
argued that the size of the dipole generating the gluon diffractive parton density should be 
small. This would also suggest a slope closer to the case of elastic J/tp photoproduction 
rather than to one expected if the soft physics dominates. Hence, in our analysis we used 



the x-dependent parameterization for B g given by Eq. (2.8). 

In addition, in order to access the sensitivity of our results to the uncertainty in the 
choice of B g , we considered the second option, when the diffractive slopes for the gluons 
and quarks are the same 

(A.12) 



B n = B, 



q > 



where B q is given by Eq. ( A.ll ). 
Model 2 



We modified the parameterizations of Eq. ( A.9Q by adding a low-/3 piece to the Pomeron 
parton distributions, to remove the rather unnatural feature that they go negative for small 
(3. As a guide, we used the factorization theorem for diffractive processes. Since small 
values of correspond to large diffractive masses, M 2 Q 2 , the diffractive cross section 



is expected to be described by formulae derived in the triple Regge limit [58] within which 
W 2 3> M 2 3> Q 2 . One assumes in this approach that multipomeron exchanges can be 
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neglected. Although there are no strong theoretical reasons to support this assumption the 
triple Regge formulae with an effective Pomeron exchange seem to work pretty well up to 
rather large energies relevant for shadowing at x > 10~ 4 . In the triple Regge approximation 
the ratio of the diffractive cross section at t = to total inclusive cross section in DIS should 
be the same as the ratio of diffractive and total cross section in real photon or hadronic 
diffraction. This may be expressed approximately (neglecting small corrections due to 
small deviations of ajp from 1) as 

d 4 a di ff/dxdQ 2 dtdMj [ \ t= o _ A 

d 2 a/dxdQi ~ Ml ' ( } 

where the coefficient A is process- independent, with a good accuracy. Indeed, data on real 
photon and pion diffractive dissociation on hydrogen give A = 0.122 ± 0.006 GeV -2 for 
the photon and A = 0.118 ± 0.006 GeV -2 for the pion [g^]. An earlier experiment on 
pion, kaon and proton diffractive dissociation on hydrogen produced similar values of A: 
A = 0.102 ± 0.013 GeV" 2 for pions at p lab = 200 GeV/c, A = 0.094 ± 0.005 GeV' 2 for 
protons at p hab = 200 GeV/c, and A = 0.097 ± 0.027 GeV" 2 for kaons at p tab = 100 GeV/c 



|6C]. In our analysis, we use A = 0.12 GeV 2 . 

Assuming this low-/? tail appears once the quark distributions become zero, using 



Eq. ( A.13j ) we arrive at the following for f3 < (3q = 0.137: 



M^o)~^ (A.14) 



and use Eq. (A. 9) for larger (5 > /3q. For the gluon channel the soft contribution is a small 



correction hence we ignore it and implement Eq. (A.9) for j3. 
Model 3 

The light-cone QCD model of Hautmann, Kunszt, and Soper |^5| leads to a set of the 
Pomeron parton densities which, upon QCD evolution in Q 2 , give a good description of 
the ZEUS data. The quark {q = u = u = d = d = s = s) and gluon parton distributions in 
the Pomeron, at the initial scale Qo = 1-5 GeV, read 



Pf q /p(0, Qo) = 0.0278/? (1 + 0.824/3 - 0.286/3^ - 2.713/3 d + 1.218/3 4 
Pfg/r(P, Qo) = 0- 987 f 1 + 0.821/3 - 1.495/3 2 + 1.569/3 3 - 0.239/3 4 ) . (A.15) 



One can see from Eq. ( |A.15| ) that, in contrast to the parameterizations of Eq. ( |A.9| ), the 
low-/3 region is treated in a sensible fashion. In this model, the Pomeron flux, /jp, is given 
by the following expression 

. , . 1 /0.0042\0.253 

fp{xip) = —( . A.16 

Xp\ Xjp J 



Note also that since this Pomeron flux is different from the one used in 17], a direct 



comparison of the overall numerical coefficients in Eq. ( |A.9j ) and Eq. ( [A. 15 ) is not possible. 
Model 4 
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The HI collaboration performed the QCD analysis of their own data [^] and produced 
an independent set of fj/jp- The computer code with this parameterization of fj/jp is 
available from fl~8[| . One should note that this parameterization is only valid in the range 
0.04 < (3 < 1. Outside of this range, the code would give parameterizations flat in j3. For 
illustrative and practical purposes, we fitted the quark and gluon parameterization to a 
simple polynomial form with a fair accuracy: 



/3/ 9 /ip(/3,Qo) 
0f a /r(P,C%) 



0.003976 - 0.02708/3 + 0.3284/3 2 - 0.7832/3 3 + 0. 



0.6105 - 1.709/? + 9.873/3 2 - 26.49/3 3 + 31.71/3 4 - 13.87/? 5 
The Pomeron flux in this model was taken in the following form 



0.3267/3 5 , 
. (A.17) 



-1.923081 4b 886 



0.01005184b 406 ) 



-2.292, 
Xjp ( 



-8.84615 + ln(xip)) . (A.18) 



B. Parameterization of the nuclear one-body density 





Po (fm" 3 ) 


c (fm) 


12 C 


0.013280 


2.2486 


32 s 


0.0049717 


3.3663 


40 Ca 


0.0039769 


3.6663 


110p d 


0.0014458 


5.308 


197 Au 


0.0008073 


6.5157 


208p b 


0.0007720 


6.6178 



Table 2: The parameters enter- 
ing the nuclear one-body density, 



p A (b,z), of Eq. (B.l) 



The nuclear one-body density pA(b, z), which appears 
in Eqs. (2.2), (2.4), (2.6), (2.14)-(2.16), (3.3), (4.1) and 
(4.7) was parameterized in a two-parameter Fermi form 



PaQ>,z) 



po 



1 + exp [(r — c)/a] 



(B.l) 



where r = y \b\ 2 + z 2 and a = 0.545 fm and the param- 
eters po and c are presented in Table 2. Also note that 
PA(b, z) is normalized as 2tt f Q °° d\b\ dz\b\pA(b, z) = 1. 
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